How robust is the entanglement entropy - area relation? 
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We revisit the problem of finding the entanglement entropy of a scalar field on a lattice by tracing 
over its degrees of freedom inside a sphere. It is known that this entropy satisfies the area law - 
entropy proportional to the area of the sphere - when the field is assumed to be in its ground state. 
We show that the area law continues to hold when the scalar field degrees of freedom are in generic 
coherent states and a class of squeezed states. However, when excited states are considered, the 
entropy scales as a lower power of the area. This suggests that for large horizons, the ground state 
entropy dominates, whereas entropy due to excited states gives power law corrections. We discuss 
possible implications of this result to black hole entropy. 



o 
o 



> 

o 



o ■ 

I 



X 



PACS numbers: 04.60.-m,04.70.-s,04.70.Dy,03.65.Ud 

Although classical black holes (BHs) have infinite en- 
tropy and zero temperature, Bekenstein - inspired by the 
area increase theorem of general relativity - proposed 
that BHs have entropy proportional to horizon area Ajj. 
This, together with Hawking's discovery that BHs radi- 
ate with the temperature Th = hc^ / (SttGM) have given 
rise to the Bekenstein-Hawking area law for BH entropy: 



Sbh — 



/ KG 

(Planck length). (1) 



The area (as opposed to volume) proportionality of 
BH entropy has been an intriguing issue for decades. At- 
tempts to understand this problem can be broadly clas- 
sified into two classes: (i) those that count fundamental 
states such as Z?-Branes and spin-networks, which are 
supposed to model BHs [3] , and (ii) those that study en- 
tanglement entropy [4, 5] and its variants such as the 
brick- wall model and Shakarov's induced gravity [6]. 

In the case of entanglement entropy, which is of interest 
in this work, it is assumed that the von Neumann entropy 



-Tr[p\n{p)\ 



(2) 



of quantum fields due to correlations between the exterior 
and interior of the BH horizon, accounts for black hole 
entropy. Such correlations imply that the state of field, 
when restricted outside the horizon, is mixed, although 
the full state may be pure [4, 5]. Although this entropy 
is ultra-violet divergent, a suitable short distance cut-off 
[0(i?p)] gives S (X Ah (it was argued in [4] that this en- 
tropy must be formally divergent). This idea gained fur- 
ther credence, when it was shown that even for Minkowski 
space-time (MST), tracing over the degrees of freedom 
inside a hypothetical sphere (of radius R), gives rise to 
the entropy of the form 0.3 {R/aY where a is the lattice 
spacing [4, 5] (it was shown in [7] that quantum fluctua- 
tions inside the sub- volume scale as its bounding area as 
well). Thus, the area-law may be a direct consequence of 
entanglement alone. 

However, a crucial assumption was made in the analy- 
ses of [4] and [5], that all the harmonic oscillators (HOs) 



- resulting from the descretization of the scalar field - 
are in their ground state (GS). Thus the natural ques- 
tion which one would ask is: How sensitive is the area 
law to the choice of the quantum state of the HOs? 

In a recent paper, the current authors had investigated 
this problem for a simpler system of two coupled oscil- 
lators, and found two interesting results [8]: (i) the en- 
tropy remains unchanged if the GS oscillator wave func- 
tions are replaced by generalized coherent states (GCS), 
and (ii) the entropy could increase substantially (as much 
as 50%) even if one of the oscillators is in its first ex- 
cited state (ES). This raises the possibility that for the 
more interesting case of iV-coupled oscillators (modeling 
a free scalar field), deviations from the area law could 
result if excited states are taken into account. We ad- 
dress this issue in this work. When the oscillators are in 
GCS and a class of squeezed states (SS), we show ana- 
lytically that entanglement entropy exactly equals that 
of the ground state, implying that the area law remains 
valid. For ESs, of the form of superpositions of a number 
of wave-functions, each of which has exactly one HO in 
the first ES, we show numerically that the entanglement 
entropy still scales as a power of the area, but that the 
power is now less than unity. The more terms there are 
in the superposition, the less is this power [9]. 

The Hamiltonian for a free scalar field (</?) is 



H = 



(3) 



where tt is the momentum conjugate of (f. Decomposing 
if and TT in terms of real spherical harmonics {Zim) i. e.. 



and discretizing on a radial lattice 



[with Tj+i 



2 



Vj = a = M ^ and L = {N+ l)a is the box size], we get, 



TV 



H, 



Im 



= -E 



•fir. 



flm,j + l 

i + 1 



(4) 



of the Hamiltonian of TV coupled HOs (up to the overall 
factor of 1/a, which does not change the entanglement 
entropy to be computed) : 



N N 



ijXj . 



(5) 



where fimjlnm,]] = 'Pimirj)[n-im.,jirj)] and [ipim,j ,^^1'^' ,y] = 
i^u'Smm'Sjj' (Note that the momenta tt;,,, j in Eq.(4) arc a with the interaction matrix elements Kij given by 
times the discretized versions of -Kim) ■ This is of the form 1, . . . , iV): 



l{l + 1) 5ij + - 5n5ji + JV - SiNSjN + + o 



A general eigenstate of the Hamiltonian (5) is given by: 
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tp{xi,. .. ,xn) ='[[Ni H^. (k^ ^ expl-^kf) ^ xj] , where TV^ = 



i 



7rV4 



(6) 



(7) 



Kd = UKU^ is a diagonal matrix {U'^U = In) with elements fc/jj, x = f/a;, Ct = U^Kl)U, such that |Q| = I^Cxfl^, 
= (xi. . . . , .T^/v), x"^ = (xi. • • • jXat) and Ui {i = \ . . . N) are indices of the Hermite polynomials {Hv). Note that the 
frequencies are ordered such that koi > koj for i > j. 

The density matrix, tracing over first n of the A'' field points, is given by: 



/n 
Yl dxi ip{xi,...,x„;ti,...,tN-n) -tp^ixi, ... ,x„;t[, ... ,t'N-n) 
i=l 
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I JJdxjexp 
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X [JiVi^f,, (4iXi) xexp 



i=l 



a;'^ • O • x' 
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^^N,H^,{ki,^), (8) 



where we have introduced the notation: t,- 



l..(iV - n), i.e. a;-' = (xi, . . . ,a:„;ti, . . . ,tAr_i) = 
{x\, . . . ,Xn\t), where t = ti,. . . , tN-n. Density matrix 
(8) yields the entanglement entropy via (2). For an ar- 
bitrary excited state (7), analytically, it is not possible 
to obtain a closed form expression for p{t;t'). When all 
the HOs arc in the GS however, i.e. when = 0, Vz, a 
closed form expression of p{t;t') and the corresponding 
entropy can be found [4, 5] . Here, we will investigate two 
non-trivial states, namely GCS and then ES. 

Before we proceed with the evaluation of S, it is impor- 
tant to compare and contrast between the two non-trivial 
states and the ground state: (i) GCS, unlike GS, is not 
an energy eigenstate of HO. However, GCS and GS are 
both minimum uncertainty states, (ii) ES, unlike GS, are 
not minimum uncertainty states. If the area law holds for 
both GCS and ES (or a superposition of ES), then this 
would indicate that it is robust and unaffected by changes 
of the chosen state. If it holds for GCS and not for ES 
(or a superposition thereof), this might signal its valid- 



ity only for minimum uncertainty states; however, if the 
reverse, or some other result holds true, a simple inter- 
pretation cannot be given and more investigation would 
have to be done. 

First, let us assume that all the HOs are in GCS, i. e.. 



^gcs{^U...,Xn) 



n 






exp 



El/2 / s2 
(xj - ai) 



(9) 

where = (ai, . . . , a^) represents the N complex GCS 
parameters . Physically, the real and imaginary parts 
of ai corresponds to the classical position (a;o) and mo- 
mentum (po) of the individual HOs respectively i. e., 

ai = xo- ipo/koi- Defining, 

X = X — U~^a , dx = dx , (10) 

it follows that 



V-Gcs (a^i • • • a:jv) 



exp 



(11) 



3 



where ipoixi . . .xn) is GS wave function. Consequently, 

/n 
Y[ dxi ipGcs {xi; t) tpacs {xr, t') 

i—1 

/n 
dxi ipo {xi; t) i/jq {xi; i') 
i=i 



ini 



where 



I^Itt^-" 
Po{i;i'), 



9.= 
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(12) 



(13) 



A, B,C are n X n, n X {N - n) and {N - n) x {N - n) 
matrices respectively, /? = \B^ A~^B , j = C — (3 and 
i = ii, . . . jtN-n = Xn,---XN- Note that B, B^ (and 
hence /3) are zero if HOs are non-interacting. Comparing 
with Eq. (11) of Ref. [5], we see that po{t,t') is precisely 
the GS density matrix. That is, the GCS density matrix 
has the same functional form, albeit in terms of the tilde 
variables. Consequently, it will have the same entropy as 
well, which is found in the following way. By a series of 
transformations, (12) can be written as: 



\n\ 



-I i N-n 



\A\n 



N-n 



where VjV'^ = = diagonal, (3 = Jjj^V(3V'^jjj^, 
WPW'^ = diagonal with elements (3i and Vi e 

V = W'^{V-fV'^)^VT. Since p^cs in Eq. (14) is product 
of the (N — n), two HO {N = 2, n = 1) density matrices, 
the total entropy is the sum of the entropies [5], i. e., 



N-n 



5= ^-ln[l-e.] 
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1-e 



1 + 



(15) 

Since the eigenfunctions and eigenvalues are the same as 
those for the ground state values, the total entropy is also 
same as that of the ground state [5]: 



1=0 



5= ^(2i+l)5, =0.3(- 



(16) 



where R = a{n + 1/2), Si is the entropy for a given 

I. Strictly speaking, the upper limit of the sum should 
have been infinity, for which it is convergent. Thus, 



Thus, one can conclude that all classical states give rise 
to the area law. This is our first result. It can be eas- 
ily shown that the results continue to hold for the class 
of SS, characterized by the same squeezing parameter r. 
The equivalents of Eqs.(9) and (10) are respectively: 



^ss(a;i,...,a;jv) 



^N 

y/r X 



r^/^ exp 



dx 



El/2 2 
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\fr dx . 



(17) 
(18) 



we choose a maximum value of I 



such that 



\S{lmax) - Silraax - 5)]/S{lma. - 5) < lO'^ (Thc nu- 
merical error in the total entropy is less than 0.1%). The 
sum in Eq.(16) is convergent. In other words, the en- 
tanglement entropy follows the area law, even for arbi- 
trarily large values of the coherent state parameters aj! 



Eqs.(ll) and (12) remain unchanged (up to irrelevant 
multiplicative factors), and so docs the entropy (16). 

Next, wc assume the state to be in a linear superpo- 
sition of N wave- functions of the form in Eq.(7), such 
that in each such wave-function there is exactly one HO 
in the first ES and the rest (iV- 1) in their GS. Although 
non-trivial, we will find that this state is relatively easy 
to handle. Using (7), the wave function for this state can 
be simplified to 



iP^{xi...xn) = 



n 



Yl ^^^^ { ^m^^ ) 



= \/2 (o^if^x) '^q{xi,. . . ,XN 



(19) 



where = (ai, . . . , cat) are the expansion coefficients 
(normalization of ijji requires a-^a = 1). 
Substituting in Eq. (8), we get 



/lb 
JJ dxi [x'^ A x'^] iPo {xf, t) {xi-, t') ,(20) 



i=l 



where A is a iV x A'' matrix and is defined as 



Al I^c ^^^^ 



and A^, Ab, Ac are nxn, nx (N — n),{N — n) x {N — n) 
matrices respectively. Comparing Eqs. (12, 20), it is 
clear that the excited state density matrix is not same 
as that of the ground state. More importantly, it is not 
possible to obtain Eq. (12) from (20) in any limit of A. 
Integrating Eq. (20), we get 



p{t,t') = po{t,t')Tr{hAA-^) 



1 - i [t^K^t + t''^k-yt') + t^Kfjt' 



where we have defined 

2Ag {A-^B) - [A-^f AaA-'B 



,(22) 



A0 = 



Tr(AAA-i) 

2A^ + [A-'] ^A^A-'B - AIA-'B - [A'^] % 
rr(A,A-i) 



Before proceeding further, we note the following: (i) 

Eq. (22) is the exact density matrix for the discretized 
scalar field with any one HO in the first excited state 
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while the rest in the GS; (ii) unhke the GS (12), the ex- 
cited state density matrix (22) contains non-exponential 
terms and hence can not be written as a product of 
{N — n), 2-coupled HO systems, as in Eq.(14). Con- 
sequently, the entropy cannot be written as a sum as in 
(15). However, for the vector outside the maximum 

of = (^-^^^^ (1, 1, . . .) , (corresponding to 3a 

limits), the Gaussian inside po{t;t') in (22) is neghgible. 
Thus, when the conditions 



> for i > j. Rewriting Eq. (25), we have: 



^max tmax ^ 1 



£1 = tmax A/3 tjnax < 1 and 62 

are satisfied, one can make the approximation: 



exp 



.(23) 



Correspondingly, the density matrix (22) takes the fol- 
lowing simple form: 



p(i,t') = 



m 



Tr{\AA-^) 



X exp (^-i [Fit + t'^yt') + t^p't'^ ,(24) 

where /?' = /3 -I- A/j, 7' = 7 + A^. Note that (24) is of 
the same form as the ground state density matrix (12), 
with matrices 7 7' in the exponent, and up 

to irrelevant normalization factors. The corresponding 
entropy will then be given by Eq.(15), with the replace- 
ments P ^ f3' and 7 — > 7' in the definition of ^i. We 
tested the validity of approximation (23) numerically for 
large values of TV (TV > 60), using MATLAB [10]. The 
error in the approximation (23) was less than 0.1% for 
— l/y/o{Q, • • • 0, 1 • • • 1) with the last o columns being 
non-zero. The corresponding entropy, computed from 
the density matrix (24) was computed numerically, for 
iV = 300, n = 100, • • • , 200 and o = 10, 20, 30, 40, 50. 

Before proceeding with the results we would like to 
mention the following: For the above ES, the expectation 
value of energy is given by: 



£1 = (Hh 



1/2 

where, fc^^ are the normal frequencies of the Hamilto- 
nian Him and £0 is the ground state energy given by 



£0 



1 



N 

i=N-o+l 



1,1/2 



(25) 



1 ^ 



1/2 

Di ■ 



(26) 



Note that we have not set the GS energy to zero and, as 
mentioned earlier, the frequencies are ordered such that 



1 /I 1 \ ^ 



„l/2 



(27) 



i=N-o+l 



As mentioned before, the kj^n are in ascending order. 
Moreover, the last few terms in the second sum in Eq.(27) 
dominate. Consequently for N = 300 and o = 10 . . . 50, 
{£1 — £0) /£o ~ 0.3 . . . 0.6 i.e. the energy of these ES are 
about 30 — 60% higher than the GS energy. Note that 
these energies are in units of 1/a, with a the ultraviolet 
cut-off (the lattice spacing). Thus, if we choose the latter 
to be of the order of the Planck length, as is reasonable in 
any theory of quantum gravity, the GS and ES energies 
(as well as the energy density of the ES) are Planckian, 
where we have ignored the resulting gravitational self- 
interactions of the system. 



ln(S ) Vsln(n + 1/2)for GS 



log(S) = 1 .9904 log(n + 0.5) - 1 .2097 



ln(S ) Vs ln(n + 1/2) for = 20 



1og(S) = 1 .8451 log(n + 0.5) - 0.4571 



in(S ) Vs ln(n + 1/2) for = 40 



1og(S) = 1 .651 log(n + 0.5) + 0.61 58 



ln(S ) Vsln(n + 1/2) for o= 10 



log(S) = 1 .8957 log(n + 0.5) - 0.7458 
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ln(S ) Vsln(n + 1/2) for o = 30 



log(S) = 1 .7695 log(n + 0.5) - 0.0337 

75I ^ ^ ^ 1 
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ln(S ) Vsln(n + 1/2) for o = 50 



log(S) = 1 .51 41 log(n + 0.5) + 1 .3869 
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FIG. 1: Logarithm of GS and ES entropies versus the radius 
of the sphere (R/a) i. e. R ^ a{n + 1/2) for iV = 300 and 
100 < n < 200. We choose the maximum value of I such that 
[S{lmax) — S{lmax — 5)]/S{lmax — 5) < 10~^. The numerical 
error in the total entropy is less than 0.1%. 



In Fig.(l), we have plotted log(S') versus 
log(i?/a). From the best-fit curves, we see that 
for o = 10,20,30,40,50, 5* = 0.4744(i?/a)0-9479, 
0.6331(i?/a)"-9223^ 0.9669(i?/a)0-8848, 1.8511(i?/a)0-8255^ 
4.002(i?/a)°'^^'^i respectively. Thus, although the 
coefficient in front increases, the power decreases with 
the number of excited states, and for large enough 
areas, the GS (or closely related GCS or SS) entropy 
is greater than the ES entropy. We would thus like 
to conjecture that if the entanglement entropy of a 
superposition of the GS and ES is computed, it would 
(at least approximately) be a sum of the GS entropy (the 
area law) and the ES state entropy that we found, in 
which case, the latter can be interpreted as (power-law) 
corrections to the area law. Such corrections can be 
contrasted with entropy corrections obtained from other 
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{2l + 1)S|Vslforn =100 (21 + 1) Vs I for n =125 
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FIG. 2: Plot of the distribution of entropy per partial wave 
[(21 + l)Si] for GS (solid-curves) and ES (dotted-curves) . To 
illustrate the difference between the GS and ES (and that all 
curves can be fitted in the same graph), we have multiplied 
the GS entropy per partial wave by a factor of 5, while the 
o = 10 and o — 30 curves have been multiplied by factors of 
6 and 2 respectively in each plots. 



higher ESs. One technical problem that we anticipate 
in this case is that the density matrix will not be ex- 
pressible as the GS matrix with shifted parameters (such 
was the case for (22)), since Hn{x) ~ x". Finally, ana- 
lytical proofs of the area law, for GCS, SS and ES, which 
do not depend on the shape of the traced out volume, 
along the lines of Ref. [14] would be illuminating. We 
hope to examine these issues elsewhere. 
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sources [11]. In Fig. (2), we have plotted the entropy 
for each partial wave, {21 + 1)5*;, versus /, for N — 300 
and various values of n. For each n, we have plotted for 
o = 10, 30, 50. It can be seen that for the GS, there is a 
maxima aX I — Q, after which [21 + 1)5*; decreases. Once 
it reaches a minimum, it starts to rise again, due to the 
large degeneracy factor {21 -\- 1). For the ES however, 
a sharp maximum occurs between / = 5 and / = 30, 
depending on the parameter o. We hope to get a better 
understanding of this phenomenon in the future. 

To summarize, in this work we have computed the en- 
tanglement entropy of scalar fields, after tracing over its 
degrees of freedom inside a hypothetical sphere of radius 
R. The oscillator modes representing these degrees of 
freedom were assumed to be in GCS, SS and a super- 
position of the first ES. In the case of GCS and SS, the 
entropy turned out to be identical to that for the GS of 
the form G.3(i?/a)^, while for the ES, the entropy goes 
as a power of the area which is less than unity. 

In the light of the above results, let us discuss the impli- 
cations of our results to the BH entropy: the Bekenstein- 
Hawking area law not only tells that the BH entropy is 
proportional to area, it also gives the precise value of the 
proportionality constant [l/(4€^)]. Our analysis, for the 
MST, suggests that the constant of proportionality as 
well as the power of the area depend on the choice of the 
state of the scalar field. This raises an immediate ques- 
tion: which states determine the Bekenstein-Hawking en- 
tropy? [12, 13]. As stated before, it appears that the GS 
is most relevant, at least for large areas, while the ES 
give rise to correction terms. 

Open problems include extending our analysis to 
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